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ABSTRACT 

In  this  report  the  application  of  the  differential  game  theory  to  missile  guidance  problem  is 
considered.  Interceptor/ target  relative  kinematics  equations  for  a  3-D  engagement  are 
derived  in  state  space  form  suitable  for  implementing  the  feedback  guidance  through 
minimisation-maximisation  of  the  performance  index.  This  performance  index  is  a 
generalisation  of  that  utilised  by  previous  researchers  in  this  field  and  includes,  in 
addition  to  the  miss-distance  term,  other  terms  involving  interceptor/ target  relative 
velocity  terms.  This  latter  inclusion  allows  the  designer  to  influence  the  engagement 
trajectories  so  as  to  aid  both  the  intercept  and  evasion  strategies.  Closed  form  expressions 
are  derived  for  the  Riccati  differential  equations  and  the  feedback  guidance  gains.  A  link 
between  the  differential  game  theory  based  guidance  and  the  optimal  guidance,  the 
proportional  navigation  (PN)  and  the  augmented  PN  (APN)  guidance  is  established. 
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Differential  Game  Theory  Application  to  Intelligent 

Missile  Guidance 

Executive  Summary 


This  report  considers  the  application  of  the  differential  game  theory  to  the  missile 
guidance  problem.  The  scenario  considered  involves  engagement  between  an  attacker 
(interceptor/ pursuer)  and  a  target  (evader),  where  the  objective  of  the  former  is  to 
execute  a  strategy  (manoeuvre)  so  as  to  achieve  intercept  with  the  target,  while  the 
objective  of  the  latter  is  to  execute  a  strategy  (manoeuvre)  so  as  to  evade  the  attacker 
and  avoid  or  at  least  delay  the  intercept.  Differential  game  approach  enables  guidance 
strategies  to  be  derived  for  both  the  attacker  and  the  target  so  that  objectives  of  both 
parties  are  satisfied. 

Interceptor/ target  relative  kinematics  equations  for  a  3-D  engagement  are  derived 
in  state  space  form  suitable  for  implementing  feedback  guidance  laws  through 
minimisation-maximisation  of  the  performance  index  incorporating  the  game  theory 
based  objectives.  This  performance  index  is  a  generalisation  of  that  utilised  by  previous 
researchers  in  this  field  and  includes,  in  addition  to  the  miss-distance  term,  other  terms 
involving  interceptor/ target  relative  velocity  terms.  This  latter  inclusion  allows  the 
designer  to  influence  the  engagement  trajectories  so  as  to  aid  both  the  intercept  and 
evasion  strategies.  Closed  form  expressions  are  derived  for  the  Riccati  differential 
equations  and  the  feedback  gains  that  allow  the  guidance  of  the  interceptor  and  the 
target  to  be  implemented.  A  link  between  the  differential  game  theory  based  guidance 
and  the  optimal  guidance,  the  proportional  navigation  (PN)  and  the  augmented  PN 
(APN)  guidance  is  established. 

The  game  theory  based  guidance  technique  proposed  in  this  report  provides  a 
useful  tool  to  study  vulnerabilities  of  existing  weapons  systems  against  current  and 
future  threats  that  may  incorporate  'intelligent'  guidance  or  for  enhancing  the 
capability  of  future  weapons  by  implementing  (game  theory  based)  intelligent 
guidance.  Further  research  is  required  in  this  area  in  order  to  evaluate  the  performance 
of  the  game  theoretic  guidance  in  realistic  combat  environments. 
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is  the  x-position  of  vehicle  i  in  fixed  axis, 
is  the  y-position  of  vehicle  i  in  fixed  axis, 
is  the  z-position  of  vehicle  i  in  fixed  axis, 
is  the  x-velocity  of  vehicle  i  in  fixed  axis, 
is  the  y-velocity  of  vehicle  i  in  fixed  axis, 
is  the  z-velocity  of  vehicle  i  in  fixed  axis, 
is  the  x-acceleration  of  vehicle  i  in  fixed  axis, 
is  the  y-acceleration  of  vehicle  i  in  fixed  axis, 
is  the  z-acceleration  of  vehicle  i  in  fixed  axis, 
is  the  x-position  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  y-position  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  z-position  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  x-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  y-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  z-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  x-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  y-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  z-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis. 

is  the  (3x1)  position  vector  of  vehicle  i  in  fixed  axis, 
is  the  (3x1)  velocity  vector  of  vehicle  i  in  fixed  axis, 
is  the  (3x1)  acceleration  vector  of  vehicle  i  in  fixed  axis, 
is  the  (3x1)  position  vector  of  vehicle  i  w.r.t  j  in  fixed  axis, 
is  the  (3x1)  velocity  vector  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  aj  -  a  j  :  is  the  (3x1)  acceleration  vector  of  vehicle  i  w.r.t  j 
w.r.t  in  fixed  axis. 

is  a  (6x1)  relative  state  vector  between  vehicle  i  w.r.t  j 


A: 

B: 

S: 

I: 


R2: 


defined  w.r.t.  the  fixed  axis. 

is  a  (6x6)  state  coefficient  matrix. 

is  a  (6x3)  control  (input)  coefficient  matrix. 

is  a  6x6  final  state  function  penalty  weighting  matrix. 

is  a  3x3  unity  matrix. 

is  a  6x6  state  function  penalty  weighting  matrix. 

is  a  3x3  pursuer's  demanded  acceleration  function  penalty 

weighting  matrix. 

is  a  3x3  evader's  demanded  acceleration  function  penalty 
weighting  matrix. 
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V  (aj1,  §2):  is  the  optimum  performance  index/ objective  function. 

H  (y  jaJjajj:  is  the  Hamiltonian 

P :  is  the  Riccati  matrix  equation  solution. 

T  =  (tf  -t):  is  the  time-to-go. 

K?,K?:  are  Interceptor  (pursuer)  state  feedback  guidance  gains. 

Ke2,K^:  are  Target  (evader)  state  feedback  guidance  gains. 

S(t)  =  n(T):  is  the  Riccati  vector  equations  solution. 
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1.  Introduction 


TACTICAL  missiles  have  been  in  use  since  WWII  and  their  guidance  systems  have 
progressively  evolved  from  those  employing  proportional  navigation  (PN)  and 
augmented  proportional  navigation  (APN)  to  those  employing  optimal  guidance  (OG)  and 
game  theoretic  guidance  (GTG).  One  reason  for  this  development  is  the  fact  that  the 
implementation  hardware  for  the  guidance  system  has  evolved  and  now  offers  greater 
flexibility  to  the  system  designer  to  implement  advanced  algorithms  for  missile 
navigation,  guidance  and  control.  Developments  in  the  area  of  IR  /RF  missile-borne  strap- 
down  seekers,  strap-down  navigation  instrumentation,  and  air-borne  processors  have 
prompted  the  guidance  engineers  to  explore  techniques  that  are  more  suited  for  the 
evolving  and  relatively  more  complex  battlefield  scenarios.  With  the  advent  of  state 
estimation  techniques  such  as  the  Kalman  Filter  and  others,  it  is  now  feasible  to  implement 
the  OG,  GTG  and  GTG  +AI  'intelligent'  guidance  on  practical  weapon  systems. 

It  is  noteworthy  that  the  PN  and  APN  are  still  being  used  in  a  number  of  modern 
missile  guidance  systems  and  give  good  performance.  PN  and  APN  performance  has  been 
studied  by  anumber  of  authors  [1-5].  It  is  also  interesting  to  note  that,  as  shown  in  [5],  both 
PN  and  APN  guidance  can  be  derived  using  the  optimum  control  theory  and  state  space 
representation  of  the  interceptor  and  target  kinematics.  Thus  we  may  regard  both  PN  and 
APN  as  special  cases  of  OG  and  GTG;  this  connection  is  further  explored  in  this  report. 
The  desire  to  reduce  weapon  life-cycle  cost,  and  at  the  same  time  extend  the  operational 
envelope  to  cope  with  complex  engagement  scenarios,  that  require  the  capability  to  adapt 
to  adversary's  'intelligent'  engagement  tactics,  it  is  necessary  to  consider  OG  and  GTG 
guidance  approaches  for  future  tactical  missiles.  Augmentation  of  these  guidance 
techniques  with  those  that  have  evolved  in  the  field  of  "intelligent  systems"  (e.g.  AI, 
Fuzzy-Neural  Nets,  Bayesian  and  others),  also  need  to  be  considered. 

Application  of  the  differential  game  theory  to  missile  guidance  has  been  considered 
by  a  number  of  authors  [6-13].  Shinar,  et.al.,  [7]  presented  an  analysis  of  a  complex  combat 
scenario  involving  two  parties  (two  aircraft),  both  equipped  with  interceptor  missiles.  The 
objective  of  each  party  was  to  shoot  down  the  opponent's  aircraft  without  their  own 
aircraft  being  intercepted  (hit)  by  their  opponent's  missile.  Such  a  scenario,  where  the 
strategy  of  each  party  is  for  its  missile  to  intercept  the  opponent's  aircraft  and  perform 
evasion  manoeuvres  of  its  own  aircraft  so  as  to  avoid  being  hit  by  opponents'  missile,  is  a 
typical  engagement  scenario  where  the  game  theoretic  approach  to  missile  guidance  can 
be  used.  Shinar  refers  to  this  situation  as  "non-cooperative  differential  game"  which  can 
also  be  classed  as  a  "game  of  a  kind".  An  encounter  between  the  two  aircraft  (blue  and 
red),  under  the  above  conditions  results  in  one  of  the  following  outcomes: 

a.  Win  for  blue-  (red  alone  is  shot  down) 

b.  Win  for  red-  (blue  alone  is  shot  down) 

c.  Mutual  kill-  (both  red  and  blue  are  shot  down) 

d.  Draw-  (both  red  and  blue  escape) 

Shinar  goes  on  to  consider  further  combat  strategies  that  can  arise  out  of  the  above 
scenario  and  suggests  the  application  of  artificial  intelligence  (AI)  augmentation  to  the 
differential  game  guidance.  This  latter  aspect  can  be  considered  from  the  perspective  of 
augmenting  the  GTG  with  a  'rule-based'  AI  -  for  switching  performance  index  weighting 
parameters  and/or  for  applying  additional  manoeuvres  to  evade  the  pursuer.  In  this 
report  our  main  focus  will  be  on  the  formulation  and  solution  of  the  GTG  problem 
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involving  two  parties  where  the  objective  of  one  party  (pursuer)  is  to  implement  a  strategy 
to  intercept,  while  the  objective  of  the  other  (evader)  is  to  implement  a  strategy  to  evade 
the  former.  Ben-Asher,  et.al.,  [8]  considered  the  application  of  the  differential  game  theory 
to  missile  guidance  and  utilised  the  optimum  quadratic  performance  index  approach  in 
order  to  derive  strategies  (in  terms  of  guidance  acceleration  commands)  for  the  pursuer 
and  the  evader  in  a  two  party  game  scenario.  The  approach  was  based  on  defining  the 
interceptor  and  target  kinematics  in  linear  state  space  form  and  the  performance  index  as  a 
quadratic  functional  of  states  and  controls  (the  so  called  LQ-optimisation  problem).  The 
above  authors  considered  engagement  kinematics  in  2-D  and  the  performance  index  was 
tailored  to  that  problem;  in  this  report  we  have  generalised  the  problem  to  3-D.  The  GTG 
problem  can  be  stated  in  a  general  form  as  follows.  Given  the  following  state  space 
(relative  kinematics)  equations: 

^i2=Ayi2+B(5i-5e2)+B(?f-??);  y12(°)=y12o  (11) 

Where: 

yi2=(x12  u12)t  :  is  the  relative  state  (position  and  velocity)  vector  between  the 

interceptor  1  and  target  2  defined  w.r.t.  the  fixed  axis. 

Xjj  =  (  Xjj  y jj  z i j ) T  :  is  the  (3x1)  position  vector  of  vehicle  i  w.r.t  j  in  fixed  axis. 

-ij  =  (uij  vij  w i j ) T  :  is  the  (3x1)  velocity  vector  of  vehicle  i  w.r.t  j  in  fixed  axis. 

(a)’ ,a2):  are  the  commanded  acceleration  vectors  respectively  of  the  interceptor  1 
(pursuer)  and  the  target  2  (evader). 

(a)1 , a 2 ):  are  the  additional  (pre-specified)  acceleration  vectors  respectively  of  the 
interceptor  and  the  target.  These  are  included  to  admit  additional  evasion  and/or  pursuit 
manoeuvres  that  the  parties  might  implement. 

The  guidance  problem  is  to  compute  the  interceptor  and  the  target  accelerations ( a j,a2j. 


such  that  the  optimum  performance  index  V*[a2,a2  j  is  given  by: 


o 


S  ,  Q  :  are  positive  semi-definite  matrices  that  define  the  penalty  weightings  on  the  states. 
Ri  ,«2  :  are  positive  definite  matrices  that  define  the  penalty  weightings  on  the  inputs. 

The  authors  [8]  solved  the  problem  (1.1)  -  (1.3)  for  a  special  case  of  2-D  with  a 
performance  index  consisting  of  miss  distance  squared  term  only.  In  this  current  report  we 
shall  consider  engagement  in  3-D  (azimuth  and  elevation  plane  engagement)  and  define  a 
performance  index  that  incorporates  the  miss  distance  term  as  well  as  additional  terms 
consisting  of  relative  (velocity)  states  that  may  allow  us  to  shape  the  engagement  trajectory 
and  more  effectively  deal  with  large  heading  errors,  unfavourable  engagement  geometries 
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and  severe  interceptor  and  target  manoeuvres.  It  must  be  pointed  out  that,  in  general,  the 
game  outcome  depends  on  which  of  the  parties  'plays  first'  [7,  8],  We  therefore  assume 
that  the  game  evolves  such  that  both  parties  apply  'strategies'  (guidance  commands)  that 
satisfies  the  'saddle  point'  solution  for  the  performance  index;  that  is  both  parties  must 
utilise  optimum  strategies,  in  which  case: 


V*(aj ,  a  2  )=  Min  Max  V  (aj1 ,  §2)=  Max  Min  V  (aj ,  a  2  ) 


(1.4) 


Both  OG  and  GTG  are  derived,  based  on  defining  an  objective  function  (performance 
index  PI)  that  is  a  function  of  system  states  and  controls  that  have  to  be  minimised  (w.r.t  to 
the  pursuer's  guidance  commands)  and  maximised  (w.r.t  the  evader's  guidance 
commands).  In  the  case  of  the  OG,  which  can  be  regarded  as  a  special  case  (a  subset)  of  the 
GTG,  it  is  assumed  that  the  target  does  not  implement  the  optimum  evasion  strategy  while 
the  interceptor  implements  the  optimum  intercept  strategy  and  does  not  perform  any 
additional  manoeuvres.  Thus,  with  a  slight  modification,  the  objective  function  (1.3)  may 
be  used  to  derive  the  OG  law.  The  optimisation  problem  for  linear  systems  that  use 
quadratic  objective  function  (the  LQ-control  problem)  was  considered  by  Sage,  A.P  [14]; 
the  solution  to  this  problem  was  obtained  through  minimisation  of  a  Hamiltonian  type 
function.  This  yields  a  linear  state  feedback  control  law  involving  'gain'  terms  derived 
from  the,  well  known,  Riccati  differential  equations  (see  section  3).  This  approach  has  been 
adopted  in  this  report  to  solve  the  GTG  and  the  OG  problems.  In  both  cases,  closed  form 
solution  of  the  Riccati  equation  and  the  resulting  feedback  gains  are  derived  as  a  function 
of  time-to-go.  While,  the  emphasis,  in  this  report  is  on  two  party  games  involving  one 
interceptor  against  one  target,  the  development  of  the  kinematics  equations  and  the 
guidance  law  derivation  is  general  enough  to  be  extended  to  multi-party  game  situation 
[7,10]. 

Intercept  and  evasion  strategies,  implemented  by  the  parties  involved,  are  based  on 
their  knowledge  of  relative  states  (i.e.  the  parties  learn  from  the  environment)  and  on 
optimisation  of  the  objective  function  (i.e.  the  decision-making  criteria).  Rule  based  AI 
methodology  can  also  be  implemented  that  will  allow  adaptively  changing  the  PI 
weightings  and/or  implementing  additional  manoeuvres  (a)1  ,§2).  In  view  of  these 
attributes,  we  shall  refer  to  the  GTG  as  'intelligent'  guidance. 

Section  2  of  this  report  presents  the  development  of  a  3-D  engagement  kinematics 
model  in  state  space  form.  Section  3  of  the  report  presents  the  formulation  of  the  objective 
function  and  the  solution  of  the  GTG  problem.  Solutions  of  the  matrix  Riccati  differential 
equations  MRDE  and  the  vector  Riccati  differential  equations  VRDE  are  considered  in 
section  4,  along  the  feedback  implementation  of  the  guidance  law.  Relationships  between 
the  OG  and  GTG  and  the  conventional  PN  and  APN  guidance  are  explored  in  sections  5 
and  6.  Section  7  contains  conclusions  resulting  from  the  material  presented  in  this  report. 


2.  Development  of  the  Engagement  Kinematics  Model 

Typical  two  vehicle  engagement  geometry  is  shown  in  Figure  2.1  (T  is  the  target  and  I 
is  the  interceptor);  this  scenario  may  be  extended  to  the  case  of  m  targets  and  n 
interceptors;  that  is  N  =  n  +  m  vehicles.  The  engagement  kinematics  model  in  this  report  is 
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derived  with  this  in  mind.  The  motion  (kinematics)  of  each  vehicle  can  be  described  by  a 
set  of  first  order  differential  equations  representing  states  of  the  vehicles  (i.e.  position, 
velocity  and  acceleration)  defined  in  fixed  (e.g.  inertial)  frame.  The  kinematics  equations 
may  be  written  as: 


d 

- X: 

dt 

d 


=  uj  ; 


—  Ui  =  av. 


dt 


(2.1) 

(2.2) 


Where: 

The  above  variables  are  functions  of  time  t . 
xj  :  is  the  x-position  of  vehicle  i  in  fixed  axis, 
y  i  :  is  the  y-position  of  vehicle  i  in  fixed  axis. 

Zj  :  is  the  z-position  of  vehicle  i  in  fixed  axis, 
u  j  :  is  the  x-velocity  of  vehicle  i  in  fixed  axis. 

V;  :  is  the  y- velocity  of  vehicle  i  in  fixed  axis, 
w  j  :  is  the  z- velocity  of  vehicle  i  in  fixed  axis. 

ax.  :  is  the  x-acceleration  of  vehicle  i  in  fixed  axis. 

ay.  :  is  the  y-acceleration  of  vehicle  i  in  fixed  axis, 

a  ..  :  is  the  z-acceleration  of  vehicle  i  in  fixed  axis. 

zi 


Note  that  a  'flat  earth'  assumption  is  made  and  that  Z-axis  is  assumed  positive  'down'. 

2.1  Relative  Engage  Kinematics  of  n  vs  m  Vehicles. 


We  now  consider  relative  state  variables  of  the  vehicles;  the  kinematics  equations  may 
be  written  as  follows: 

d  d  d 

- X::  =11::;  - y::  =  Vjj;  - Z  ::  =WH 

dt  IJ  1J  dt 3  J  ,J  dt  ,J  IJ 
d  d  d 

—  uj:  =  a...  =  ax.  -a... ;  —  v.-.-  =  av..  =  av.  -av. ;  — w=  =  a7..  =  a7.  -a7. 

dt  J  M  1  J  dt  J  y*J  Yl  yJ  dt  'J  1  J 


(2.3) 

(2.4) 


Where: 


XU 

y>j 

Zij 


u 


IJ 


vij 


W: 


a... 

x'j 

av.. 

y'J 


=  X;  -  Xj  :  is  the  x-position  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  yj  -yj  :  is  the  y-position  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  zj  -Zj  :  is  the  z-position  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  uj  -Uj  :  is  the  x-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  Vj  -  vj  :  is  the  y-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  w;  -  w  j  :  is  the  z-velocity  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  ax.  -ax.  :  is  the  x-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  ay.  -a  v.  :  is  the  y-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis. 

=  az.  -aZj  :  is  the  z-acceleration  of  vehicle  i  w.r.t  j  in  fixed  axis. 
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T 


Line  of  Sight 


I:  Interceptor 
T:  target 


Figure  2.1- Interceptor/Target  Engagement  Geometry 


In  the  context  of  the  engagement  between  an  interceptor  against  a  target  (2-party 
engagement),  we  may  regard  suffix  i  to  represent  the  interceptor  and  suffix  j  to  represent 
the  target.  In  the  derivation  of  the  optimal  guidance  law  it  will  be  useful  to  represent  the 
above  equations  in  vector/ matrix  notation. 

2.1.1  Vector/Matrix  Representation 

We  can  write  equations  (2.1),  (2.2)  as: 

—  x.  =u.;  —  u.=a.  (2-5) 

dt  -1  -1  dt  -1 

Similarly  we  can  write  the  relative  kinematics  equations  (2.3),  (2.4)  as: 


d 


d 


(2.6) 


a 


-J 


Where: 

x .  =  (  X;  ys  z  j )  1  :  is  the  (3x1)  position  vector  of  vehicle  i  in  fixed  axis. 

Uj  =  (  Uj  Vj  Wj )  1  :  is  the  (3x1)  velocity  vector  of  vehicle  i  in  fixed  axis. 

=  ( ax.  ay.  az. ) 1  :  is  the  (3x1)  acceleration  vector  of  vehicle  i  in  fixed  axis. 

Xy  =  (  Xjj  y  y  z i j ) T  :  is  the  (3x1)  position  vector  of  vehicle  i  w.r.t  j  in  fixed  axis. 

Ujj  =  ( u jj  v jj  wy ) T  :  is  the  (3x1)  velocity  vector  of  vehicle  i  w.r.t  j  in  fixed  axis. 

a-  =  (ax„  ay„  az„  )T  =  as  -  a  j  :  is  the  (3x1)  acceleration  vector  of  vehicle  i  w.r.t  j  in  fixed 

axis. 
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Equation  (2.6)  may  be  combined  together  to  give  us: 


V 

'o 

1 

V 

+ 

"o' 

a  j  - 

"o' 

ij  _ 

0 

0 

_-ij_ 

I 

I 

(2.7) 


This  can  be  written  as: 

—  y  =Ay  +Ba;-Ba; 
dt  -i.i  ->j  -1  -J 


(2.8) 


Where: 

y  =  ( X;j  Ujj)T  :  is  a  (6x1)  relative  state  vector  between  vehicle  iand  vehicle  j  defined 
w.r.t.  the  fixed  axis. 

'o  r 


A  = 


B  = 


0  0 
0 
I 


is  a  (6x6)  state  coefficient  matrix. 

:  is  a  (6x3)  control  (input)  coefficient  matrix.,  I  is  the  identity  matrix. 


Remarks 

In  this  report  the  guidance  algorithm  is  derived  based  on  linear  engagement 
kinematics  defined  in  the  fixed  (e.g.  inertial  axis).  There  are  thus  implied  assumptions  in 
the  derivation:  (a)  that  the  body  attitude  does  not  vary  significantly  during  the 
engagement,  and  (b)  that  the  guidance  commands  are  implemented  in  the  fixed  axis. 
However,  in  most  engagements  the  body  attitude  may  vary  significantly  during  the 
engagement  and  more  importantly  the  guidance  commands  are  applied  in  the  body  axis 
and  require  the  guidance  commands  (generated  in  fixed  axis)  to  be  transformed  to  body 
axis.  For  testing  and  performance  of  the  guidance  strategies  derived  in  this  report,  a  non¬ 
linear  engagement  kinematics  model  should  be  used  for  the  simulation  model.  Such  a 
simulation  model  was  derived  by  the  author  and  is  given  in  [11],  which  will  be  referred  to 
as  the  Integrated  Navigation,  Guidance  and  Control  Test  Bed  (INGC-Test  Bed).  The 
guidance  commands  should  be  applied  in  the  vehicle  body  axis  (through  appropriate 
transformation)  which  accounts  for  changes  in  body  attitude;  also  the  autopilot  dynamics 
and  effects  such  as  the  gravity  and  vehicle  lateral  acceleration-limits  should  be  included  in 
the  simulation  model.  Several  authors  [5]  have  included  autopilot  lags  in  the  guidance  law 
derivations  and  the  material  presented  in  this  report  may  be  extended  to  that  case. 


3.  Optimum  Interceptor/Target  Guidance  for  a  2-Party 

Game  Scenario 


We  shall  first  of  all  address  the  problem  of  a  target  that  is  being  engaged  by  an 
interceptor  and  wishes  to  implement  a  guidance  strategy  to  avoid  intercept;  while  the 
interceptor  endeavours  to  implement  a  strategy  to  effect  an  intercept.  The  respective 
strategies  are  implemented  via  the  application  of  guidance  acceleration  commands.  A 
scenario  of  this  type  will  be  referred  to  as  a  2-party  game. 
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The  interceptor  i  =  1  utilises  its  guidance  input  command  a2  to  effect  the  intercept  of 
the  target,  this  is  specified  by  ;  in  addition  it  may,  if  required,  perform  an  additional 
pre-specified  manoeuvre  a  2  .  The  total  interceptor  acceleration  in  this  case  may  be  written 
as: 

aj  =  aj1  +aj  (3-1) 

The  target (j  =  2)  on  the  other  hand,  applies  an  evasive  manoeuvre  a2 ,  and  an 
additional  pre-specified  (disturbance)  manoeuvre  a2  >  this  latter  component  could  be,  for 
example,  a  random  manoeuvre  or  a  manoeuvre  of  a  periodical  wave  form.  With  these 
manoeuvres,  the  total  evader  acceleration  is  of  the  form  (see  equations  (2-8)): 

^2  ~~2  +  -2  (3-2) 

The  modified  form  of  the  kinematics  equation  (2.8),  including  the  above  manoeuvres 
may  be  written  as: 

ftln=Ali2+B(^+^)-B(^+^)  <3-3) 

The  object  of  the  2-party  game  is  to  compute  the  evasion  and  pursuit  guidance 
(strategies)  commands  a2  ,  aj1  that  satisfy  specified  criteria.  For  this  report  we  shall 
consider  the  application  of  the  differential  game  and  the  optimum  control  principles  to 
derive  both  the  evasion  and  pursuit  guidance  commands  (strategies). 

3.1  Construction  of  the  Differential  Game  Objective  Function 

In  the  context  of  missile  guidance,  differential  game  ensues  if  the  interceptor's 
objective  is  to  utilise  a  strategy  (i.e.  execute  guidance  commands)  that  results  in  an 
intercept  with  the  target,  while  the  objective  of  the  target  is  to  utilise  a  strategy  (i.e. 
executes  guidance  command)  so  as  to  delay  or  totally  avoid  intercept  and  eventually 
escape  from  the  interceptor.  Here,  intercept  implies  that  the  minimum  separation  between 
the  target  and  the  interceptor  is  small  or  preferably  zero.  In  formulating  the  DG  guidance 
problem,  the  following  assumptions  are  made: 

(a)  Both  parties  have  all  the  necessary  information  of  the  states/ relative  states,  with 
respect  to  each  other,  to  enable  the  parties  to  implement  the  necessary  guidance 
laws.  Countermeasures  designed  to  conceal  states  of  parties  involved  are  not 
considered. 

(b)  If  it  is  assumed  that  a  seeker/ target  tracker  is  used  to  construct  system  states  from 
seeker  information  (e.g.  utilising  a  Kalman  Filter).  However,  for  the  purpose  of  our 
current  considerations  it  is  assumed  the  system  states  (or  state  estimates)  are  exact. 
Later  studies  will  include  a  state  estimator. 

(c)  The  maximum  and  the  minimum  accelerations  (both  lateral  and  longitudinal) 
commanded  and  consequently  achievable  by  the  vehicles  involved  in  the  game  are 
limited.  For  our  current  derivation  constraints  on  the  accelerations  are  considered  to 
be  "soft";  that  is,  the  change  in  acceleration  and  acceleration  rate  is  gradual  in  the 
neighbourhood  of  maximum/ minimum  values.  This  type  of  constraint  can  be 
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implemented  in  the  objective  function  through  the  use  of  'penalty  weightings' 
associated  with  the  demanded  acceleration,  and  leads  to  a  relatively  easier  solution 
and  implementation  of  the  guidance  law. 

(d)  In  the  derivation  of  the  guidance  laws,  autopilot  lag  is  ignored.  However,  this  may 
be  included  in  the  actual  simulation  studies  in  order  to  assess  the  guidance 
performance  including  the  autopilot. 

Under  the  above  assumptions  we  can  proceed  to  construct  the  objective  function 
(performance  index)  that  the  parties  to  the  game  will  need  to  minimise  or  maximise  in 
order  to  derive  their  respective  strategies.  Note  that  the  interceptor's  objective  is  to 
minimise  the  relative  separation  (miss  distance),  while  that  of  the  target  is  to  maximise  the 
miss  distance.  The  objective  function,  therefore,  must  include  interceptor /target  states  that 
represent  the  miss  distance  as  well  as  other  states  that  influence  this.  In  the  past,  most 
authors  (8,  10)  have  used  miss  distance  (which  is  a  function  of  relative  position)  and 
demanded  accelerations  terms  to  construct  the  objective  function.  In  this  report  we  shall, 
generalise  the  objective  function  by  including  terms  in  relative  position  and  relative 
velocity  (states)  as  well  as  demanded  accelerations  (input).  Thus  the  objective  function 
may  be  written  as: 


o 


Where: 

v(-)  :  is  the  objective  functional. 

(x12 1  x  12  j:  is  the  square  of  the  relative  range  (separation)  between  the  interceptor  and  the 
target.  (x12  1  u12  ):  is  the  projection  of  the  relative  velocity  on  to  the  relative  range. 

(uj2  1  u12  ):  is  the  square  of  the  relative  velocity. 

:is  the  square  of  the  interceptor  (pursuer's)  acceleration. 


:  is  the  square  of  the  target  (evader's)  acceleration. 


v  J 

(s!,S2,S3):are  penalty  weighting  on  the  final  values  of  state  functions. 
( q , ,  q  2 1  <|  3 ) :  are  penalty  weighting  on  state  functions. 

( q ,  r2 ) :  are  penalty  weighting  on  demanded  accelerations. 


By  varying  the  relative  values  of  the  penalty  weightings  (s^q^q),  constraints  on 
system  states  and  control  can  be  implemented.  Note  that  the  objective  function  contains: 

(a)  Interceptor/ target  relative  position  term:  (x121x12);  at  the  final  time  t  =  tf  it  is  the 
miss  distance  squared. 

(b)  Interceptor/ target  relative  position  and  velocity  terms:  (x12Tu12)  and  (u l2 1  u]2] 
these  terms  represent  engagement  trajectory  shaping  terms. 
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(c) 


Interceptor/ target  demanded  acceleration  terms: 


pT  p 
“i  »i 


and 


these  terms 


v  y  v  j 

allow  soft  constraints  on  controls  (demanded  accelerations)  to  be  implemented. 

(d)  Relative  values  of  penalty  weightings:  (  s1,s2,s3),(  q1;q2,q3)  and  (  rhr2)  determine 
soft  constraint  boundaries  on  states  and  control  variables. 

The  objective  function  v(..)  as  given  in  equation  (3.4)  can  be  minimised  w.r.t  ajj*  in 


order  to  derive  the  intercept  guidance  commands,  and  maximised  w.r.t  ,  by  virtue  of 
the  negative  sign  associated  with  it,  in  order  to  derive  the  evasion  guidance  commands  for 
the  target.  It  will  be  convenient  to  write  the  objective  function  as: 


(  p  A  1 

2 

\  y12»»i*»2j=2 

Zi2 

s 

(3.5) 


Where: 

2sjl  s2I 
s2I  2s3I 
matrix. 

2qii  q2! 

q2l  2q3I 

Rl  =riI:  is  a  3x3  pursuer's  demanded  acceleration  function  penalty  weighting  matrix. 

K  2  =  r2I  :is  a  3x3  evader's  demanded  acceleration  function  penalty  weighting  matrix. 

ll  ll2  ta 
|  a  ||  =  a  A  a  . 


Q  = 


■  is  a  6x6  final  state  function  penalty  weighting  matrix.  I :  is  a  3x3  unity 


:  is  a  6x6  state  function  penalty  weighting  matrix. 


In  this  report  we  will  consider  the  objective  function  with  Q  =  0 .  The  game  theoretic 
guidance  problem  can  be  stated  as  that  of  minimising  the  objective  function  V |yj2,aj  ,§2) 

w.r.t  aj*  and  maximising  it  w.r.t  a2  ,  that  is: 


Min  Max  Vl 
„P 


(-l2’-l ’-^)= 


Min  Max 
HP  -2 


In 


tf  , 

-,*5  (I" 


2 

Ri 


2 

R2 


dt 


(3.6) 


Remarks 

For  a  minimum  and  a  maximum  of  the  objective  function  to  exist,  it  is  a  requirement 
that  the  matrix  S  be  at  least  positive  semi-definite,  and  matrices  k  (  and  R2  be  positive 

definite.  That  is,  the  determinants:  |s|  >  0  and  |R  j  |  >  0 ,  |R2 1  >  0 .  These  conditions  imply  that 
(see  Appendix  A.l): 

sl5s2>0,and  (4sjs3 -s22)  >  0  ,  s2  can  be  positive  or  negative  (3.7) 

3.2  Solution  of  the  Differential  Game  Guidance  Problem 


In  this  section  we  present  the  solution  to  the  problem  of  optimisation  of  the  objective 
function  (3.6)  subject  to  the  condition  that  the  kinematics  equation  (3.3)  holds.  In  doing  so 
we  shall  follow  the  technique  described  in  [14].  This  involves  the  construction  of  the 
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Hamiltonian,  which  is  minimised  w.r.t  and  maximised  w.r.t  a2  .  The  Hamiltonian  may 
be  written  as: 

af-  si  )=  i  {(rfTRi*r)  -  (a5TR2*l)}  -  •• 

+  -(A^12  +  B(-1  +-l)-B(-2  +  -2)) 

Sufficient  conditions  for  Min  Max H  (yi2,aj,a2)  are  given  by: 


5?  -2 


(a) 

5  H( 

5a  j  ' 

,yi2'5i,a^)=0  and 

(b) 

8  H 

(-12’-1 ’-^)=  ~~  — 

(c) 

The  terminal  condition  for  : 

5a® 


(3.8) 


H2' 


The  corresponding  necessary  conditions  are 

(a)  -^-TH(yi2,aP,a*)>0  and  -^-H  <  0 

5af  5a® 


Now,  applying  the  sufficiency  conditions  (a)  to  the  Hamiltonian  (3.8),  we  get: 


-^-n(yu,3:pl,ae2)=R13P1+BTX  =0 

5a  j 

(3.9) 

^H(y  aP,a®)=-R2a®2-BTk  =  0 

5a  2 

(3.10) 

Equations  (3.9),  (3.10)  give  us: 

a\  =-R1“1BTk;  and  a®=-R2_1BTk 

(3.11) 

It  is  left  to  the  reader  to  verify  the  necessary  conditions.  We  are  interested  in 
constructing  guidance  commands  that  are  functions  of  system  relative  kinematics  states; 
we  shall  therefore  assume  that  h  is  of  the  form: 

k  =  py12+|  (3.12) 

Where: 

P :  is  a  6x6  matrix,  which  will  be  later  shown  to  be  the  solution  of  the  matrix  Riccati 
differential  equation  (MRDE),  and 

is  a  6x1  vector,  which  will  be  later  shown  to  be  the  solution  of  the  vector  Riccati 
differential  equation  (VRDE). 

Thus: 

§1  =-Rr1BTPyi2  -Rf'B1^;  and  a®  =  -R^^P  y^  -R^B^  (3.13) 
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We  shall  write  these  equations  as: 

=  ~K1  y12  and  ae2  =  -K*  yJ2  (3.14) 

Where: 

Kj  =  R!_1BTP  :is  the  state  feedback  gain  for  the  pursuer;  and  Kj  =R1-1B1  :  is  the 
disturbance  input  gain  of  the  pursuer. 

Ke2  =R21BtP  :  is  the  state  feedback  gain  for  the  evader;  and  K2  =R2  'b1  :  is  the 
disturbance  input  gain  of  the  evader. 

Now  applying  the  sufficiency  condition  (b),  we  get,  using  equation  (3.8): 

-^H(yi2,aP,ae2)=ATL  =  4  (3.15) 

—12 

Where: 

X  =  P  +P  +4 
-  -12  -12  5 

Substituting  for  X  ,  X  ,  and  y  ,  and  for  ajj*  and  a2 

simplification  (detailed  derivation  is  given  in  Appendix  A.2) 
equation  (3.15)  leads  to: 

{p  +  PA  +  AtP-Pb(r,_1-R2'1)btP  }  y  ••• 

=  {_4_[A-B(Rf1-R2-1)BTp]T-PB(aJ-a^)} 

Since  the  solution  of  equation  (3.15)  must  hold  for  all  y  ,  a  solution  of  this  equation 


can  be  written  as: 

P  +  PA  +  ATP-Pb(r1“1  -R2_1)  btp  =  o 

(3.17) 

And: 

-^-[a-bIr!-1  -R2_1)  BtP  JT-PB(af  -a2)=0 

(3.18) 

With  terminal  conditions  P(tf)=S  ,  and  tf)  =  0  .  In  the  sequel,  equation  (3.17)  will 
be  referred  to  as  the  MRDE  and  equation  (3.18)  will  be  referred  to  as  the  VRDE. 

4.  Solution  of  the  Riccati  Differential  Equations 

4.1  Solution  of  the  Matrix  Riccati  Differential  Equations  (MRDE) 

In  this  section  we  shall  consider  deriving  the  solution  of  the  MRDE  (3.18).  We  write 
these  equations  as: 

P  +  PA  +  AtP-PBRBtP  =  0  (4.1) 

Where: 

r  =  (r1-1-R2“1) 


,  and  after  algebraic 
it  can  be  shown  that 

(3.16) 
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Note  that:  for  R1=2r1I,  R2=2r2I, 
r2  >  r3  .  In  fact  if  r2  »  rj  ^  R  =  2r3  I . 


R  =  2 


r  r  r  ^ 

rlr2 
r2  -U 


I ;  and  for 


Rto  be  positive  definite 


The  approach  adopted  to  solve  the  MRDE  involves  an  inverse  matrix  technique,  where 
the  solution  an  inverse  matrix  version  of  the  MRDE  is  first  obtained  and  then  by  re¬ 
inverting  the  resulting  solution,  the  Riccati  matrix  [P]  is  obtained.  This  approach  is  given 
in  Appendix  A.3.  Equation  (4.1)  is  solved  in  Appendix-3  in  some  detail  to  obtain  an 
expression  for  E  and  then  re-inverted  to  obtain  expressions  for  elements  of  P .  For  the 
particular  case  considered  in  Section.3,  (where:  sn  =s22  =  s33  =  2s3;  s14  =s25  =s36  =s2; 


S44  =s55  =s66  =2!*3  and  rli  =r22  =r33  =2r,T  =  tf  -  t)  we  get: 


24r 

4rsj  + 

4*1*3  -s22 

T 

48r 2  +  48rs3  T  +  24rs2  T2 

+  16rs j  T3  + 

4SiS3  -S22 

H 

■t*. 

1 _ 1 

(4.2) 


P  44  -  P  55  -  P  66  ~ 


P 14  -  P25  -  P26 


8r  [l2rs3  +12rs2T  +  12rsiT2  +(4sjs 

3-22 

T3] 

48r 2  +  48rs3  T  +  24rs2  T2  +16rsj  T3  + 

12r  [4rs2  +8rs1T  +  (4s1s3  -s2 

,4s1s3 

2)  T2] 

-u 

1  T4 

48r 2  +  48rs3  T  +  24rs2  T2  +16rs!  T3  + 

4sis3 

-22! 

x4 

(4.3) 


(4.4) 


4.1.1  Pursuer/Evader  State  Feedback  Guidance  Gains 


The  feedback  gain  matrix  for  the  interceptor  is  given  by: 


Kf  =—  BTP  =  — 

2rj  2rj 

P  14 

0 

0 

P44 

0 

0  ' 

0 

P25 

0 

0 

P55 

0 

(4.5) 

0 

0 

P  36 

0 

0 

P  66  _ 

The  feedback  gain  matrix  for  the  target 

is  given  by: 

K  j  =—  BTP  =  — 
2r2  2r2 

P  14 

0 

0 

P44 

0 

0  ' 

0 

P25 

0 

0 

P55 

0 

(4.6) 

0 

0 

P  36 

0 

0 

P  66  _ 

4.1.1.1  A  Special  Case 

Another  case  of 

interest  is 

when: 

■11  =s 

22  =  s33  =  2sU 

s14  =s25  =s36  =s2 

S44  =  s55  =  s66  =  2s3  =  0  and  r31  =  r22  =  r33  =  2r  .  This  is  equivalent  to  setting  the  weightings 
on  the  velocity  terms  in  the  performance  index  to  zero;  in  which  case  the  latter  is  only  a 
function  of  the  miss  distance  squared  and  terms  in  interceptor  and  target  accelerations. 
For  this  case: 

6rs.  ..  _ 

P11  =  P22  =  P33  =  T  r  (4.7) 

[3r  +  Sl  T3] 


12 


UNCLASSIFIED 


UNCLASSIFIED 


DSTO-TR-2838 


6rsjT 

P 14  =  P  25  =P36  =T - 

3r +  Sj  T3 


P44  -P55  -P66  - 


6rsjT2 
3r  +  Sj  T3 


P44 

P14I 

Pll 

P55 

=  T 

P25 

_  rp’  ^ 

P22 

_P66_ 

_P36  J 

.P33. 

(4.8) 

(4.9) 


(4.10) 


Substituting  for  pjj  from  (4.8)  -  (4.9),  the  feedback  gain  matrix  for  the  pursuer  in  this  case 
is  given  by: 


1  0  0  T  0 

K**=—  BtP=  ..  3rSlT — =-  0  1  0  0  T 

2ri  rJsr  +  SjT3]  0  0  10  0 


0 

0 

T 


(4.11) 


The  feedback  gain  matrix  for  the  evader  is  given  by: 


1  0  0  T  0  0 

„  1  T  3rsiT 

Kj  = - 8^  = — f - I - =-0  1  0  0  T  0 

2r2  r2[3r  +  SiT3J  0  0  1  0  0  T 


(4.12) 


4.2  Solution  of  the  Vector  Riccati  Differential  Equations  (VRDE) 

The  VRDE  given  in  equation  (3.18)  may  be  written  as: 

4  =  -  [a-BR_1BTP  ]T§-PB  (ad-ad)  (4.13) 


Writing:  ^  =  [  t,2  ^3  %4  ^5  5e] T ;  equation  (4.13)  (in  its  decomposed  form)  may 


be  written  as  (see  Appendix  A.4)): 

si = c, -p  14  (•;,  -*;2)  (4.14) 

(«;,-•«)  <415> 

(<— J2)  <416> 

i4=-Sl+^-S4-P44(<,  -<2)  (4.17) 

is=-S2+— S5-P!5(=“  -*!,)  (4-18) 

r22  v  yi  y2/ 

i6  =  -^+m^-P66^dz ,-a"  )  (4.19) 

r33  '  1  2 


2  to  obtain  analytical  solutions  to  equations  (4.14)  - 

2:  (ad  ,ad  ,ad  ,  (ad  ,ad  ,ad  )  and  (ad  ,ad  ,ad  ), 
V  u  xi  V  >i  >i  '\>  \  zj’  z}’  zj/' 


Unfortunately  it  is  not  easily  possible  to  obtain  analytical  solutions  to  equations  (4.14) 
(4.19),  except  for  special  cases  where: 
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i  =  1,2  are  constants.  This  case  will  be  considered  later  on  in  this  section.  In  general, 
however,  equations  (4.14)  -  (4.19)  have  to  be  solved  backwards  in  time.  For  this  purpose 
we  make  the  substitutions: 


LetT  =  tf-t,*  dT  =  - dt ;  Jj (t)  =  $(tf -T)=  t,  (t);  aj.(t)=aj.(tf -T)=ad(T);i  =  l,2  y  =  x,y,z. 
Hence  the  above  equations  (4.14)  -  (4.19)  may  be  written  as: 

dTii 


dT 

dq2 

dT 

dr|3 

dT 

dri4 


Pl4  |  d 

- h4-Pl4 

rll  '  1 


P  25 
r22 
P  36 
r33 


95  —  P  25 

96  -P36 


P44 

—  —  T|i  H - T|4  —  P44 

dT  1  rn  4  44 


-H2  +  — ^5 
r22 


d95 
dT 

^  =  -93  + 
dT  r33 


(ad  -ad  ) 

(ad  -ad  ) 

(ad  -ad  ) 

(ad  -ad  ) 

\  xi  X2  / 

(ad  -  ad  ) 
V  yi  y 2  / 


P55 


-adJ 


(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 


These  equations  satisfy  the  boundary  condition  that  9  (d)=  (tf  )  =  0 ,  and  must  be 

solved  backwards  in  time,  that  is,  T  ->  0  .  In  the  absence  of  an  explicit  closed  form  solution 
to  equations  (4.20)  -  (4.25)  we  shall  write  the  general  solution  as: 


(T)-V  T,pjj,(axi  ax2),  (ayi  ayz),  (azj  az2) 


4.2.1  Analytical  Solution  of  the  VRDE 


(4.26) 


Analytical  solution  of  the  VRDE  is  possible  for  the  case  when:su  =s22  =s33  =2s2; 
S14  =s25  =s36  =«2  =0;*44  =  s5s  =  s66  =  2s3  =  0  and  ru  =r22  =r33  =2r;  ax.,ay.,az.  are 

constants;  equations.  For  this  case  (see  Appendix-A.4),  the  solution  of  equations  (4.20)  - 
(4.25)  gives  us: 


9  l 


9  2 


9  3 


9  4 


9  5  = 


3rs,T 


3r  + 

SjT3 

[  3rs 

rp  2 

1  1 

1 

+ 

s  jT3 

1 

(/} 

rp  2 

l1 

+ 

U 

_ i 

SjT3 

i 

-s 

C/5 

rp  3 

l1 

1 

+ 

SjT3 

3rs 

rp  3 

1  1 

3r  +  SjT' 


(ad  -  ad  ) 
(ad  -  ad  ) 
(ad  -  ad  ) 
(ad  -  ad  ) 

( an  -aJ2) 


(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 
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4  6 


3rsjT' 


(ad  -  ad  ) 

\  zi  Z2  / 


(4.31) 


3r  +  SjT' 

Noting  that:  r|  (T)  =  \  (tf  -  T),  the  feedback  for  the  disturbance  term  may  be  written  as: 

(ad  -ad  ) 

V  XI  x2/ 


KdT1  =  Rr1BTr1  =  -1- 
Zri 


4  4 

1 

4  5 

2rj 

4  6. 

3rsxT3 
3r  +  2s!T3 


(ad  -ad  ) 

V  yi  J'2  > 

(ad  -ad  ) 

\  zi  Z2  / 


Kdr1  =  R2-1BT11  =  -1- 
-  2r2 


4  4 

4  5 
4  6 


1 

2r2 


3rs!T3 
3r  +  2SJT3 


(ad  -ad  ) 

(ad  -ad  ) 

V  yi  >2  > 

(ad  -ad  ) 

V  zi  z2  / 


(4.32) 


(4.33) 


4.3  Mechanisation  of  the  Game  Theoretic  Guidance 


Using  the  expressions  for  the  optimum  guidance  equations  (4.32)  and  (4.33),  and  those 
for  the  feed-back  gains  given  in  equations  (4.32)  and  (4.33),  we  get,  for  the  general  case: 

S11  =  s22  =  s33  =  2sl  r  s14  =  s25  =  s36  =  s2  )  s44  =  s55  =  s66  =  2s3  and  rll  =  r22  =  r33  =  2r  : 

Pl4x12  +  P44u12  +  44 

P25Y12  +P55X12  +hs  |  (4-34) 

LP36z12  +P66w12  +rl6. 

P  14x12  +  P44u12  +  44 

P25Y12  +  P55v12  +  45  |  (4.35) 

I  P 36z  12  +  P66w12  +  46 


a?  =  -K?  y  -  K?  ri  = - — 

-1  1  -12  1 J  2rj 


a2=-Ke2y 


2£_12 


■K2h=-^- 

2r2 


This  is  depicted  as  a  network  block  diagram  Figure  4-1.  For  the  case  when: 

S11  =  s22  =s33  =2si;  s14  =s25  =s36  =s2  =°is44  =s55  =s66  =2s3  =0  and  rll  =  r22  =  r33  =  2r ) 

ax.  ,ay.,az.  constants,  we  get: 


^i=-Kh12-Kin= 


-  SrsjT 


Xii  +T1117  +  T 


(ad  -  ad  ) 
\  xj  X2  / 

- r - 1 —  y Tvi?  (cid  -  ad  ) 

f  V  yi  y2' 

2r1[3r  +  2s1  T3J  .  . 

(ad  -ad  ) 

V  Zi  Z2  / 


Z17  +  Tw  1 7  +  T 


(4.36) 


UNCLASSIFIED 


15 


DSTO-TR-2838 


UNCLASSIFIED 


—2 


-Ke2yi2-K^  =  - 


-  3rsjT 


2r2[3r  +  2s j  T3] 


x12  + 
yi2  + 


Tu12+T2(a^-<2) 

Tvn+T^-ajJ 

(ad  -  ad  ) 

\  Z1  z2  / 


z  12  +  Tw12  +  T 


(4.37) 


5.  Extension  to  Optimum  Guidance 


The  development  of  the  optimum  guidance  not  involving  target  evasion  manoeuvre 
can  proceed  directly  from  the  game  theory  development  presented  in  earlier  sections.  For 
this  case,  the  objective  function  takes  the  following  form: 


(5.1) 


S  = 


-  is  a  6x6  final  state  function  penalty  weighting  matrix.  I :  is  a  3x3  unity 


Where: 

'2s  j  I  0 

0  °_ 

matrix. 

Ri=ril:  is  a  3x3  pursuer's  demanded  acceleration  function  penalty  weighting  matrix. 


Note  that  in  this  case,  the  first  term  in  the  objective  function  is  only  the  miss  distance 
and  that  the  target  game  theory  based  evasion  manoeuvre  is  not  present.  Similarly  in  the 
kinematics  equations  pre-specified  target  acceleration  a2  is  present  whereas  ad  is  zero. 
The  kinematics  equation  for  the  engagement,  equation  (3.3)  may  be  written  as: 


—  y  =Ay  +BaP-Ba? 
dt  -12  -12  -1  —2 


(5.2) 


Following  the  approach  presented  earlier  it  follows  that  the  optimum  guidance  law  for 
the  interceptor  in  this  case  is  given  by: 

= -Ri_1BTP  yj2 -Rr‘BT^  = -Kj  y12  -Kd  \  (5.3) 


Or  in  terms  of  time-to-go  (5.3)  has  the  form 

aP  =  -Rf^P  -Ri_1BTrj  =  -Kj  y^  -K?  q 


(5.4) 


Where  the  MRDE  and  VRDE  are  the  same  as  those  previously  derived  in  section  4.  That  is: 


P  +  PA  +  ATP-PBR1"1BTP  =  0 

(5.5) 

And: 

-^-[a-BRj^B1?  ]T+PBad  =0 

(5.6) 
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As  far  as  the  solution  of  these  equations  (5.5),  (5.6)  is  concerned,  these  are  identical  to 
those  derived  earlier  with  r  replaced  by  rj  and  s2=s3=0  It  can  be  easily  verified  that,  for 
this  case: 


6riSi 

Pll  -P22  -P33  -  r  1 

(5.7) 

[3jr  +  Sj  T3 J 

6rjSjT 

Pl4  -P25  -P26  -  r  I 

(5.8) 

[3r1+SlT3] 

OrjSjT2 

P44  -P55  -P66  -  r  1 

(5.9) 

3rx  +sj  T3] 


The  feedback  gain  matrix  for  the  interceptor  in  this  case  is  given  by: 

1  0  0  T  0  O' 

0  1  0  0  T  0 

0  0  1  0  0  T 

And  the  disturbance  term  for  constant  target  manoeuvre  a!)  is  given  by: 


d 

1 

B4 

B? 

.Be. 

1 

3s!T3 

a 

x2 

ad 

y2 

ad 

L  z2  J 

2rj 

~~2 

3rj  +2sxT3 

RP  =—  BtP  = 
•  1  - 


3si  T 


2r 


[3r1+SlT3] 


(5.10) 


(5.11) 


Remarks 

It  is  interesting  to  note  that  if  rx  -»  0,  then  (5.11),  (5.12)  become: 


K 


—  BtP  =  — 

2r,  T2 


1  0  0  T 
0  10  0 
0  0  10 


0 

T 

0 


0 

0 

T 


K?  r|  =  Ri_1BTr|  =  — — 
-  2r, 


B4 

a*2 

3 

d 

Bs 

~~2 

ay2 

.Be. 

ad 

L  z2  J 

(5.12) 


(5.13) 


This  relationship  is  the  'link'  between  the  optimal  guidance  and  the  conventional 
guidance  such  as  PN  and  APN,  and  will  be  further  elaborated  in  the  next  section. 


6.  Relationship  with  the  Proportional  Navigation  (PN) 
and  the  Augmented  PN  Guidance 

In  order  to  establish  a  connection  between  the  optimal  guidance  and  the  PN  and  APN 
we  shall  assume  that  the  engagement  trajectory  is  such  that  the  azimuth  and  elevation 
sightline  angles  y  21, 0  21  (Figure  A-l)  remain  small  during  engagement,  that  is,  the 
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trajectory  remains  close  to  collision  close  geometry.  For  this  condition  it  follows  that  the 
interceptor/ target  relative  velocity  is  pointed  approximately  along  the  sight  line  and  is 
approximately  equal  to  the  closing  velocity  Vc .  Using  equations  (5.10)-  (5.13),  the  state 

feedback  guidance  acceleration  aj'  we  shall  write  as: 


j.p 


K5’y12-K?!i=-3Vc 


VCT 


VCT 


x12  + 


-yi2  + 


Vc  — yzi2  + 


VCT 


VCT 


VCT 


“12 


v12 


"12 


3 

+  — 
2 


a 


x2 

ad 

yi 

d 

t 

*2 


a 


(6.1) 


It  is  shown  in  Appendix  A5,  equations  (A5-3)-(A5-5)  that: 
f  x  f  \ 


¥21 


Vl2  .  yi2 


vVcT  ycT2 


;  02i  - 


^12 


"12 


VCT  Vt2 


;  and  —^xn  +  -^-u12  =  0 


VcT 


VCT 


Thus  equation  (6.1)  reduces  to: 


aP=3Vc 


'  0  " 

ad  ’ 

x2 

3 

d 

¥21 

+  — 

2 

a 

yi 

02i 

ad 

L  z2  J 

(6.2) 


This  is  the  well-known  augmented  proportional  guidance  (APN)  law  when  the  target 
manoeuvre  is  constant.  The  navigation  gain  is  3VC  associated  with  the  sightline  rate,  and 

3 

—  associated  with  target  acceleration.  When  there  is  no  target  manoeuvre  we  get  the 
proportional  navigation  (PN)  guidance. 


7.  Conclusions 


This  report  was  focused  on  the  application  of  GTG  where  the  interceptor's 
(pursuer's)  objective  was  to  hit  the  target  while  the  target's  (evader's)  strategy  was  to 
avoid  intercept.  Since  the  parties  involved  in  the  engagement  derive  their  respective 
strategies  based  on  the  assessment  of  the  current  relative  states  (i.e.  the  environment)  and 
optimisation  of  the  objective  function  (i.e.  the  decision  making  criteria),  we  shall  term  this 
type  of  guidance  as  ‘intelligent’  guidance.  It  was  further  shown  that  OG  was  a  special  case 
of  the  GTG.  Both  these  guidance  techniques  follow  similar  procedures  for  deriving  the 
state  feedback  guidance  algorithm;  the  objective  function  (performance  index),  however, 
in  the  two  cases  were  different.  A  3-D  interceptor/ target  engagement  was  considered  and 
the  various  guidance  laws  were  derived  that  will  allow  design  engineers  the  flexibility  to 
choose  guidance  gains  (by  selecting  appropriate  performance  index  weights)  to  meet  their 
specific  engagement  objectives.  For  a  number  of  important  cases  closed  form  expressions 
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have  been  obtained  for  the  feedback  gains.  A  relationship  between  the  GTG,  the  OG  and 
the  classical  PN  and  APN  has  also  been  demonstrated. 

The  game  theory  based  guidance  technique  proposed  in  this  report  may  provide  a 
useful  tool  to  study  vulnerabilities  of  existing  weapons  systems  against  current  and  future 
threat  systems  that  may  incorporate  'intelligent'  guidance  or  for  enhancing  the  capability 
of  future  weapons  by  implementing  (game  theory  based)  intelligent  guidance.  Further 
research  is  required  in  this  area  in  order  to  evaluate  the  performance  of  the  game  theoretic 
guidance  in  realistic  combat  environments.  Further  research  is  required  to  test  the 
guidance  technique  proposed  for:  imperfect  measurement  of  the  guidance  seeker,  system 
lags,  and  robustness  to  variations  in  parameters  and  noise. 
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T 

Figure  5-1.  Network  Block  diagram  of  the  guidance  mechanisation 
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Appendix  A 

A.l  Verifying  the  Positive  Semi-Definiteness  of  Matrix  S 

Now,  the  determinant  of  the  block  matrix  S  may  be  written  as: 

ISu  S12| 


|s|  = 


>12  °22 


=  s 


11 


s22-Si2Ts1r1s12 


(Al-1) 


For  S  to  be  positive  semi-definite,  we  must  have: 


lsl  = 


2sjl  s2I 

1 

S22 

1 

f  s221 

II 

bJ 

2s3I — ?-I 

II 

Ki 

CZ5 

2s3  — — 

s2I  2s3I 

2s  i 

2sj 

^  J 

>0 


(Al-2) 


That  is:  s3  >0,  and 


2s ; 


**2 

2s 


2\ 


>0  or  4s3s3  >  s22 


v  *  y 

Note  that:  s2  can  be  either  positive  or  negative. 


(Al-3) 


A.2  Derivation  of  Riccati  Differential  Equations 


Substituting  for  a)'  and  from  equations  (3.13)  into  equation  (3.3)  gives  us: 

.-•if) 


=  Ay„  -B  (Rr‘BTPy12  +  Kf'll'c 


dt  —12  J-\2 

i 


+  b  ( R2_1BTPy  12  +  R2_1bt^ — §2 ) 


(A2-1) 


Equation  (3.15)  may  be  written  as: 

aT(  py,2 + ^)=  -  i = -p  y  ,2  - p  i,,  -  i 


Substituting  for  y  from  equation  (A2-1)  gives  us: 


pyi2-piAyi2-B 


(py12+|)— Py,2-f{ 

( 


(R,-|BTPy„+RI-1BT 


+  B(R21BTPyi2+R21BT| 


—  12 

-i"))  -5 


l-sf) 


ATPy1!+AT5  =  -Py|1-PAyi2+PB(R,-'BTPy„+R1“'BT5-a 


-12  -12 


-12 


-if)’ 


PB(R21BTPyi2+R2-1BT^-a^)  -§ 


(A2-2) 


Re-arranging  the  terms,  equation  (A2-2 

[p  +  PA  +  ATP  -  PB  (r,-1  -  R2_1)  BTP 
=  -  4 -  [  aT  -  PB  (Rf1  -  R2_1)bT]  %■ ~  PB  (af  -  a*) 


can  be  written  as: 


—12 


(A2-3) 


A  solution  of  equation  (A2-3)  is  obtained  if  the  LHS  and  the  RHS  are  both  equal  to  zero 
that  is: 
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[p  +  PA  +  ATP-Pb(r1“1  -R2_1)bTp]  y  =0 

And 

-4-[at-Pb(r1-1-R2“1)bt]  ^-PB(a^-a^)=0 


Since  the  above  equations  hold  for  all  y  ,  a4  and  a  ^  ;  we  get: 

P  +  PA  +  AtP-Pb(r1“1  -R2_1)  BtP  =  0  (A2-4) 

-4-[at-Pb(r1“1-R2'1)bT]  ^-PB  (a^  -aj  )=  0  (A2-5) 

These  differential  equations  satisfy  the  boundary  conditions  P  ( tf )  =  S  and  £,  ( t ,  )  =  0 .  In 

the  sequel,  equation  (A2-4)  will  be  referred  to  as  the  Matrix  Riccati  Differential  Equation 
(MRDE)  and  equation  (A2-5)  will  be  referred  to  as  the  Vector  Riccati  Differential  Equation 
(VRDE).  For  convenience  we  shall  write:  R-1  =  (Ri_1  -R2  *)■ 

A.3  Solving  the  Matrix  Riccati  Differential  Equation 

Let  us  write  the  P  matrix  as: 

P  =  E_1;-»  PE  =  I;-»  PE  +  PE  =  0;->  P  =  -E_1EE_1  (A3-1) 


Substituting  for  P  in  equation  (A2-4,  we  obtain  the  inverse-MRDE  for  E  as: 

E  =  AE  +  EA  T  -  BR  -1 B  T  (A3-2) 


We  shall  solve  for  E ,  the  inverse  Riccati  matrix  first,  and  then  invert  this  to  obtain  the 
Riccati  solution  P  .  Because  both  P  and  E  matrices  are  symmetric;  we  may  write: 


ell 

®12 

e13 

e14 

®15 

®16 

e12 

®22 

e23 

®24 

e25 

®26 

®13 

e23 

®33 

e34 

e35 

®36 

e14 

®24 

e34 

®44 

®45 

e46 

®15 

®25 

®35 

®45 

®55 

®56 

®16 

e26 

®36 

e46 

®56 

®66 

(A3-3) 


The  terminal  condition  for  matrix  E(»f)  is  given  by:  E(tf  )  =  S  1 .  Let  us,  therefore  define: 


Sn  S12 
Si2  s22 


Sjj  0  0  s14  0  0 

0  s22  0  0  s25  0 

0  0  s  33  0  0  s36 

s14  0  0  s44  0  0 

0  *25  0  0  s55  0 

0  0  s36  0  0  s66 


(A3-4) 
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Where: 


[S»]= 

Sjj  0  0 

0  s22  0 

;  [si2]= 

Sy  0  0 

0  s25  0 

'  [  s  22  ]  - 

S  44  0  0 

0  s55  0 

0  0  S33 

0  0  s36_ 

0  0  s66_ 

For  the  particular  case  considered  in  section  3.1,  sn  =  s22  =  s33  =  2s2 ;  s14  =  s25  =  s36  =  s2  , 
and  s44  =  s55  =  s66  =  2s3 .  We  shall  first  consider  the  general  matrix  S  as  shown  in  equation 
and  later  substitute  the  particular  values.  We  write  the  inverse  of  the  inverse  matrix  S-1  as: 


S1  =  T  = 


ln 


‘12 


1 12 


*22 


(A3-5) 


Note  that,  since  S  is  symmetric  then  so  is  T  =  S  1  and  for  T  to  be  the  inverse  of  S,  we 
must  have: 


'S11 

S12' 

'T11 

T12 

I 

o' 

S12 

S22  . 

T.2 

H 

0 

I 

-(SiiTh+S12T12) 

(SllTl2  +  Si2T22  ) 

I  0 

_(^  12 T]i  +  S22T12  ) 

(S12T12  +S22T22  )_ 

0  1 

(SnTn+S12T12)  =  I 

(SiiTi2+S12T22)  =  0 

(Sl2Tll+S22T12)  =  0 

(S12T12  +S22T22)  =  I 

Equation  (A3-8)  4 

Tn  =_s12  1s22t12 


Substituting  equations  (A3-10)  into  equation  (A3-6)  gives  us: 
T12  =  [S12  _S11S12  1  ^ 22 ]~ 

Substituting  equation  (A3-11)  into  equation  (A3-10)  gives  us: 

Tll=-S12  *S22 [S12  _S11S12  1s22^ 

Also,  equation  (A3- 7)  gives  us: 

T22=-S12  1Si2T12 


Equations  (A3-11)  and  (A3-13)  give  us: 
T22=-Si2  1S11[s12 -SxjS^  ^22]^ 

It  follows  from  (A3-4)  that: 


(A3-6) 

(A3-7) 

(A3-8) 

(A3-9) 


(A3-10) 


(A3-11) 


(A3-12) 


(A3-13) 


(A3-14) 
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1 

0 

0 

1 

0 

0 

1 

0 

0 

*11 

*14 

*44 

0 

1 

0 

;  s12-1  = 

0 

1 

0 

>  S22  1  = 

0 

1 

0 

*22 

1 

*25 

1 

*55 

1 

0 

0 

0 

0 

0 

0 

*33  _ 

*36  _ 

*66  _ 

Using  expressions  for  Sjj  and  Sjj  1 ,  equations  (A3-11),  (A3-12)  and  (A3-14)  give: 
s44 


Tn  = 


(  *11*44  * 14  2  ) 

0 


0 

*55 


Tl2  = 


*14 


(*11*44  * 14  2 ) 

0 


(*22*55  * 25  2 ) 

0 


0 

~  *25 

(*22s55  * 25  2 ) 

0 


0 

0 

_ *66 _ 

(*33*66  —  *36  2  ) 

0 

0 


-  s 


36 


l22 


'll 


(  *11*44  * 14  2 ) 

0 


*22 


(*22*55  * 25  2 ) 

0 


(*33*66  —  *36  2  ) 

0 

0 

_ ^33 _ 

(*33*66  —  *36  2 ) 


(A3-15) 


(A3-16) 


(A3-17) 


S-1  =  T  = 


matrix  can  now  be  constructed  using  expressions  for  T]],T]2,T22. 

S66=2s3,  then 


T11  T12 
_T12  T22 

Note  that,  when  sn  =  s22  =  s33  =  2sj;  s14  =s25  =s36  =s2,  ands44 
using  equations  (A3-15)  -  (A3-17)  it  can  be  shown  that:: 

2s3I  -s2I 

s2I  2s2I 


*55 


,-l 


1 


(4*l*3  — *2  2  ) 


(A3-18) 


A.3.1.  Solution  of  the  Inverse  Matrix  Riccati  Differential  Equation 

In  this  section  we  decompose  equation  (A3-2)  in  its  elemental  form  to  facilitate  the 
solution  of  the  inverse  MRDE  for  E . 
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Now: 


1 

o 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

o' 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

* 

> 

H 

II 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

And: 


0  0  o' 

0  0  0 

ODD 

0  0  0  1  0  0 

o 

o 

B  = 

u  u  u 

inn 

;*bt  = 

0  0  0  0  1  0 

;  and  we  write  R  = 

0  rn  0 

A  U  U 

0  1  0 

0  0  0  0  0  1 

o 

o 

0  0  1 

(A3-19) 


It  follows  that  from  the  above  that: 


e14 

e24 

e34 

e44 

e45 

e46 

e14 

e15 

e16 

0 

0 

o' 

e15 

e25 

e35 

e45 

e55 

e56 

e24 

e25 

e26 

0 

0 

0 

e16 

e26 

e36 

e46 

e56 

e66 

;  EAt  = 

e34 

e35 

e36 

0 

0 

0 

0 

0 

0 

0 

0 

0 

e44 

e45 

e46 

0 

0 

0 

0 

0 

0 

0 

0 

0 

e45 

e55 

e56 

0 

0 

0 

0 

0 

0 

0 

0 

0 

_e46 

e56 

®66 

0 

0 

0 

Also: 


'o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1bt  = 

0 

0 

0 

1/ 

0 

0 

An 

0 

0 

0 

0 

1/ 

0 

A22 

0 

0 

0 

0 

0 

A 

Hence  the  RHS  of  equation  (A3-2)  may  be  written  as: 

AE  +  E  AT-BR'¥ 


2ei4 

(e15  +e24) 

(e16  +  e34) 

e44 

e45 

e46 

(e15  +e24) 

2e25 

(e26  +  e3s) 

e45 

e55 

e56 

(e16  +e34) 

(e26  +e35) 

2e36 

e46 

e56 

e66 

e44 

e45 

e46 

_  1/ 

An 

0 

0 

e45 

e55 

e56 

0 

-x, 

0 

e46 

e56 

^66 

0 

0 

_  1/ 
A33 

(A3-20) 


(A3-21) 


(A3-22) 
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Since  E  is  a  symmetric  matrix,  we  need  to  consider  only  the  elements  of  the  upper 
triangular  matrix.  Thus,  using  equation  (A3-22),  differential  equation  (A3-2)  may  be 
written  (in  its  elemental  form)  as: 

ejj  =  2e14;  e12  =  (e15  +e24);  e13  =(e16  +e34);  ^4=644;  e15  =e45;  e16  =e46; 

®22  =  2e25i  ®23  =  (e26  +e35)i  ®24  =  e45i  ®25  =  e55i  ® 26  =  e56' 

®33  =  2e36i  ®34  =  ®46 '  ® 35  =  e56i  ®36  =  e 66  ^ 

®44  =_/rjj  '  ^45  =  *L  ®46  =  0 


f  ®56  -  o; 


®66 =-  y 


The  terminal  conditions  are  given  by  E(tf )  =  S  1  may  be  written  as: 


ell  ( 

*f  )  = 

s44 

\  '  e22  (  tf  )-  / 

X  /  ®  33  (  L 

(  2 

V  S11S44  _s14 

j  \  s  22  ^55  —  s 

25  'j 

e14  ( 

h)=  j 

-S14 

'  e2s(tf  )—  / 

e36(tf  ) 

(  2\ 

(slls44  -  s14  / 

\  s  22  s  55  _  s  25 

2) 

®44  ( 

b)  =  ■ 

S11 

'  e55(tf  )-  / 

7T;  e  66  (  ^  f  ) 

(  2\ 

1 

(slls44~s14  / 

Vs  22  s  55  _s25 

2) 

(  s  33  s  66  _  s  36 
-s36 

t  A 

( s  33  s  66  -  s  36  / 

_ yu _ 

[  s  33  s  66  _  s  36  2  ) 


Integrating  the  above  differential  equations  with  terminal  conditions  and  writing 
!=(«,-«)  ,  time-to-go,  we  get  expressions  for  e  jj ;  this  are  given  in  Table  A-I  below. 


Table  A-l 
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Table  A-l  ( contd ) 


f  I  I 

e44 (tf  )=  e44 (t) -  Xli  I d(T  =  e 44 (t) -  f  ;  *  e44 (t)  =  e44 (tf  )  +  ^ 


e44  (t)  = 


*11 


+  CD;^e44(t)= 

2 )  rn 


( 


S11  s44  ~  s14 


,( 


1 


rns ii  +\sns44  _s 


rlllslls44  _s14 


1  1 

e3e(tf  )  =  e36(0+ J  e66(o)dor  =  e36(t)+|^  e66(tf  )+X33  (tf  -<*)  da 

t  t 

*  e36(t)  =  e36(tf)-e66(tf)(tf  -t)-^tf 


2r 


33 


■^e36(t)-_ 

^e3e(t)=- 


s36 


s33 


(tf  -t)- 


(tf  -t) 


(s33s66  _s362)  (s33s66  _s362) 

2r33s36  +  2r33s33T  +  (s33s66  ~  s362)  T2 


2r 


33 


2r33  (*33*66  “  s362) 


tf 

e35(.,)  =  e35(.)+j 

t 

e56 (a) da  =  e35 (t);  ^  e35 (t)  =  e35 (tf  )  =  0 

tf 

e34(tf  )=  e34(t)+  J 
t 

e46 (a) da  =  e 34  (t);  ^  e34  (t)  =  e 34  (t f )  =  0 

1 

e33(tf  )  =  e33(t)+2 

t 

tf 

j*e36(a)  da 

e33(tf  )-e33(t)  +  2J[e36(tf  )_e66(tf  )(tf  _t)_/4r33  (tf  _t^ 


da- 


e33(t)-e33(tf  )_2e36(tf  X*f  _t)+e66(tf)(tf  -t)2  +  (*f  -t)3; 


33 


e33(t)  = 


s66 


■  +  - 


2s36(tf-t)  s33(tf-t)2  (t 


■  + 


[s33s66  _s362)  (s33s66  _s362)  (s33s66_s362 


) 


+  — 


^e33(t)  = 


3r33s66  +  6r33s36T  +  3r33s33T2  +  ( 

* 33 s 66  _  s362  , 

1  T3 

1 

-5 

u* 

^ * 33 s 66  _  s362 

)  J 

UNCLASSIFIED 


UNCLASSIFIED 


Table  A-l  (contd) 


DSTO-TR-2838 


lt 

e  26  (t  f  )  =  e  26  (* )  + 1  e  56  (^)  =  e  26  (0 ;  ^  e  26  (t )  =  e  26  (t  t  )  =  0 


1  1 

e2s(tf  )  =  e2s(0+ J  e5s(CJ)dCT  =  e25(t)+ )+y^22  (*f  _CJ) 


da; 


e25(t)-e25(tf  )_e55(tf  X*f  _t) 


(tf-O2 


■^e25(t); 


s25 


2r22 

s22(tf-t)  (tf_t)2 


2)  (-  -  "  2 


is22s55_s25  )  \s22s55_s25 


2r 


22 


*e25(t)  = 


2r22s25  +2r22s22T  +  l 

ls22s55 

“S252, 

1  T2 

1 

-! 

,s22s55 

1 

‘1 

e24 (tf  )  =  e24  (t)+  J 
1 

• 

e 45 (°W  =  e 24(0;  ■*e24(t)=e24(tf)=0 

tj 

e23  (*f  )  =  e23  (t)+  J 
( 

• 

[e 26 (°')  +  e 35 (CT)] dcj  =  e23 (t) ;  ■*  X23(t)=X23(tt)=0 

e22  (tf  )  =  e22  (0+  J  2e25  (CT)  da 

t 

‘f 

e22  (tf  )  =  e22  (*)+  ^  2|^e25  (tf  )—  e5s(tf  )  (*f  _CT)_/2r22  (tf  _CT) 


da 


e22 (*)  —  e22  (tf  )—  2e25  (tf  )(tf  _t)+e55(tf)(tf  _t)2  +  “  “ 

*  e22  (0 


3r 


22 


2s25(tf_t)  ,  *22  (tf  —  t)  2  ,  (tf  —  t)  ' 

>  /  \  •  /  \  ' 


*  e22(t)  = 


(s22s55  _s252)  (s22s55  _  s252)  (s22s55  _  s252) 

) 


3r22s55  +6r22s25T  +  3r22s22T2  +_(f22s55  —  s  25  2  J  T3 


3r22  \  s22s55  _s25 


ei6(tf)  =  ei6(t)/  *  e16(t)=e16(tf)=0 
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Table  A-l  (contd) 


e15  (tf  )  —  e15  (*)  ^  ^  e15  (f)  —  e15  (ff  )  —  0 


el4(tf )  =  +  J  e44^)  dcr  =  e14(t)+ 1  [  e44(tf )+  Yr  (tf  -  a)]  da 


^  e14(tf)- ei4(t)+e44(tf)(tf -t)- 


(tf-Q' 


2r, 


11 


*  e  14  (0 = 


s14 


s  11  (ff —  (tf -02 


(*11*44  s  14 2 )  (*11*44  _*142)  11 


*  e14  (t)=  - 


2rll*14  +  2rll*llT  +  1 

*11*44 

S  0 
_  *14 

I  T2 

2rn  | 

*11*44 

s  2' 

“  *14 

1 

t 

f 

e13  (tf  )=£l3(tf  )+  J 
1 

[e16 (CT)+  e34  (a)]  ttCT  =  e13  (tf  );  e13  (t)=  e13  (tf  )=  0 

t 

t 

e12  (tf  )  =  eJ2  (tf  )+  J 
1 

f 

• 

[e15 (CT)+  e24  (a)]  ttCT  =  e12  (tf  );  e12  (t)  =  e12  (t f  )  =  0 

i 

4 

eii(tf)  =  en(0+J,2e14(a)da 
t 

tf 

e  1 1  (t f  )  =  £n(t)+  j" 2^e14(tf )- e44(tf  )(tf  -a)-^r^  (tf  -<0” 
t 

ell(t)=eii(tf)-2e14(tf )(tf  - 1)+  e44(tf )(tf  -t)‘  +^- — — 


da 


ll 


JV  „  ff\  *44  2*14 (tf  —  t)  Sll (tf  —  t) 2  (tf-l 

"7  rr  7  A  (  IT 

\ *11*44  —  *14  /  \ *11*44  —  *14  )  \*11*44  —  *14  ) 


*  ell(0  = 


3rns44  +6riiSi4T+3riiSiiT  +  Sjjs44— s^4  T 


3rll(*ll*44  — *14" 
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In  view  of  table  A-l,  we  may  write: 

ejj  0  0  C|4  0  0 

0  e22  0  0  e25  0 

E_  0  0  e33  0  0  e36 

e14  0  0  e44  0  0 

0  e2s  0  0  e55  0 

0  0  e36  0  0  e66 


(A3-23) 


Inversion  of  the  E  matrix  to  obtain  P  can  proceed  in  the  same  way  as,  shown  earlier 
for  S .  Since  P  is  a  symmetric  matrix,  we  write  it  as: 


E_1 


P12 

P22 


(A3-24) 


Where: 


ell 

0 

0  ' 

e14 

0 

0  ' 

e44 

0 

0  " 

En  = 

0 

e22 

0 

'  E,2  - 

0 

e25 

0 

'  e22  - 

0 

e55 

0 

0 

0 

e33  _ 

0 

0 

1 

V© 

m 

a> 

0 

0 

1 

VO 

VO 

a> 

And: 


1 

- 

i 

- 

1 

- 

0 

0 

0 

0 

0 

0 

ell 

e14 

e44 

0 

1 

0 

;  e12_1  = 

0 

1 

0 

;  s22  1  = 

0 

1 

0 

e22 

e25 

e55 

0 

0 

1 

0 

0 

1 

0 

0 

1 

e33  . 

e36  . 

e66  . 

For  P  to  be  the  inverse  of  E ,  we  must  have: 


En 

E12 

Pll 

P12  " 

I 

o' 

_E12 

E22_ 

_P12 

^0 

w 

1 _ 

0 

I 

(A3-25) 


Using  expressions  forEjj  and  E;j  1 ,  it  follows  from  (A3-25)  that: 

e44 


li 


(elle44  e  14  2  ) 
0 


0 

e55 


(e22e55  —  e25  2  ) 

0 


0 

0 

^66 


(e33e66  e 36  2  ) 


(A3-26) 
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P12  - 


_e14 


(ene44  -e142) 
0 


P22  = 


ell 


( e  1 1 e 44  _el42) 
0 


0 

_e25 


(e22e55  e 25  2  ) 

0 


0 

_ L22 _ 

(e22e55  —  e25  2  ) 

0 


0 

0 

'  e36 


(e33e66  —  e36  2  ) 

0 

0 

_ L33 _ 

(e33e66  e 36  2  ) 


Where: 


Pll 

0 

0  ' 

Pi4 

0 

0  ' 

p44 

0 

0  ' 

Pll  = 

0 

P22 

0 

;  P12  - 

0 

P25 

0 

to 

to 

II 

0 

Pss 

0 

0 

0 

P  33  _ 

0 

0 

P  36  _ 

0 

0 

P  66  . 

Substituting  for  ejj,  from  Table  A-I  we  can  now  derive  expressions  for  p 
(A3-26)  -  (A3-28).  A  detailed  derivation  is  shown  below: 

Now, 


P11  = 


nil 


e44 


■li 


i  elle44  _e14 


Where: 


nn  -  e44(0 
dn  =( 


«-l 


rl  1 s  1 1  + 

S11S44  -S14 

2 

j 

T 

rll 

sns44  -s142 

ene44  _ei4 


3rnS44  +  6rjjSj4T  +  +(^11844 

-s142)t3 

3rn(sns44  -  si42) 

2rlls14  +  2rllsllT  +  ( s  1 1 s 44  -  *142)t2 

2 

riisn  +(sns44  ~si42) 
rn(sns44  _si42) 


2ri  1  ( s  1 1 s 44  -  si42) 


(A3-27) 


(A3-28) 
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du  = 


3rii2sus44  +3r11s44(r11S44  -s142)T  +  6r112s11s14T  +  6r11s14(r11s44  -s142) 


+  3rn2Sn2T2  +4rnS||  (rns44  -  S|42 


)t3+( 


rns44  -Sn2)2  T4 


2„  2  ,  .  2„  2^2 


3rll2(silS44  _s142) 
(slls44  s  14 ~  ) 


1 4rn  *14  +  4rn  s41  T  +^S11S44 -s14' j  T  +8r14  sns14T 
1  +  4rlls14 (slls44  _s142)t2  +4rllsll(  slls44  _s142)t3 


4riii(sns44  -s142) 


dn  = 


l2rn"(slls44  —  s142  )+  l2rlls44  (  slls44  _  s142  )T  +  l2rlls14(rlls44  _s142)t2 


< 

+  4i'nS  11 1 

(rilS44  _s142) 

|t3+I 

(sll 

S44  _s142) 

I2 

1  T 

4 

12rn2| 

S11S44  -s 

2' 

14 

I2 

Pll 


III 

Ll 


121-n 

rn*n + 1 

r  ^ 

S11S44-S14 

1  T 

] 

12IJJ2  +  12rns44  T  +  12rns14  T2  +  4rjjSjj  T" 

+  ( 

*11*44  _  *142) 

|  T4 

By  proceeding  in  the  same  manner  as  above  it  can  be  shown  that: 

P22 


12r22 

r22*22  +  1 

r 

.*22*55  “  *252  , 

It 

12r222  +  12r22*55  T  +  12r22*2s  T 

+  4r22s22  T3 

♦( 

.*22*55  _  *252) 

1  T4 

P33 


12r33 

r33*33  + 

r  ' 

,  *33*66  “  *362 

T 

l2r332  +  l2r33*66  T  +  121*33*36  T 

2  +  4r33s33  T3 

♦  ( 

*33*66  ~s362, 

T4 

Now: 


P44 


11 44  _ _ ^11 _ 

1,44  .elle44  ~  e142) 


Where: 


n44  = 


3i*hS44  +  6r41s14T  + . 

7 

^rllsnT  + 

[» 

US44  -*i42 

rp3 

3rnl 

SHS44  —  *14 

1  J 

d44  - 


12ri12(s11s44  -Si42)+12r11s44(s11s44  -Si42)  T  + 12r11si4(r11s44  -Si42)  T2 


+  4r11s11(r11s44  -s142)t3  +(sns44  -s142) 


2  I^t4 


slls44-s14  )  T 


12rii2(sijS44  -S142) 


(A3-29) 

(A3-30) 

(A3-31) 
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P  44 


4rn 

3rns44  +  6rj  jS  j4T  +  3rj  jS  jjT  2  + 

s  1 1 s 44  “s14  . 

T  3 

12 Fjj  +  12rjjS44  T  4-  12rjjSj4  T2  + 

T3  + 

,*11*44  - 

S142 

)  T4 

Similarly,  it  can  be  shown  that: 

P55  = 


4r22  [3r22s55  +  ^r22s25T  +  3r22S22T2  +  S22S55  ~S252)  T3 


12r222  +  12r22S55  T  +  12r22s25  T"  +  4r22s22  T3  +  (s22s55  -s252)  T4 


P66  = 


4r33 

3r33s66  +  ^r33s36T  +  3r33S33T2  + 

*33*66  _s362, 

t3 

l2r332  +  12r33s66  T  +  12r33s36  T2  +  4r33s55 

T3  + 

,*33*66  “ 

‘s362 

)  T4 

Further: 


P 14 


Where: 


*14 


■14 


-ei4 


dl4  (elle44  _  e142  ) 


2rll*14  +  2rll*llT  + 

*11*44 

s  2) 

—  *14 

|  T2 

2rnl 

Jl2rn  2  (  slls44 

(*11*44 

-*142, 

-*,42) 

)+12rn 

1 

*44 ( *11*44  ‘ 

+  4rll*lllrll*44  _* 


14 


’44  —  *14 


'14 


')  T  +  12rlls14(rlls44  s  14 2  ) 


l2rn"(*ll*44  _si42) 


6rn 

2rll*14  +  2rll*llT  + 

,*11*44  _  *14 

2 

rj-i  2 

12 rjj  +  12rjjS44  T  +  12rjjSj4  +  4rjjSjj  + 

,*11*44 

-»,42) 

T4 

Pl4  “ 


Similarly,  it  can  be  shown  that: 

P25 


P  36 


6r22 

2r22*25  +  2r22*22T  +  (*22*55  “  *252)  T2] 

12r222  +  12r22s55  T  +  12r22s25  T2  +  4r22s22  T3  +  (s22s55  - s252)  T4 

6r33 

2r33*36  +  2r33*33T  +  (*33*66  _  s 36 2 )  T2 

42r33  2  +  l2r33* 

66  T  +  12r33s36  T2  +  4r33s33  T3  +  (s33s66  -  s362)  T4 

A.4  Solution  of  the  Vector  Riccati  Deferential  Equation 

(A2.5),  tha 

5-pb(.;-«') 


Let  us  now  consider  equation  (A2.5),  that  is: 

[  .  T  lx 


A-BR_1BtP 
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(A3-36) 

(A3-37) 
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Now: 


a-br_1btp 


0  0  0 


'  P 14 


rll 

0  0  0  0 

0  0  0  0 

-P44 


1  0  0 


'll 


0  10  0 


0  0  1  0 


~  P  25 
r22 

0 

0 

~  P  55 
r22 
0 


0 

0 

—  P  36 
r33 

0 

0 

—  P  66 

r33 


;  pb  = 


Pl4 

0 

0 

P  44 

0 

0 


0 

P25 

0 

0 

P55 

0 


0 

0 

P  36 

0 

0 

P  66 


(A4-2) 


Writing:  \  = 

^3 

^4 

^5  ^e] T  >  equation  (A4-1)  (in 

its  decomposed  form)  may 

written  as: 

4i 

=  —^4  P 14  ( 
rll  X 

ad  ) 
x2 ' 

(A4-3) 

ii 

=  ^5-pd 
r22 

(ad  - 
^  yi 

) 

(A4-4) 

43 

=  ^<U-P3„( 
r33 

ad  - 
•  zi 

ad  ) 

z2  / 

(A4-5) 

44 

=  -5,+£4L?4 

I'll 

-P44I 

(ad 

1  xj 

-<) 

(A4-6) 

4S 

r22 

—  P  55 

(ad 

V  yi 

(A4-7) 

46 

=  -63+^56 
r33 

—  P  66 

(ad 

V  ZJ 

-  ad  ) 

z2  > 

(A4-8) 

Unfortunately  it  is  not  easily  possible  to  obtain  analytical  solutions  to  equations  (A4-4) 
-(A4-9),  except  for  special  cases  where  (a  d,  a  d,  a  d],  (ad.,ad.,ad  j  and  (ad,ad,ad],  i  =  1,2 


are  constants.  This  case  will  be  considered  later  on  in  this  appendix.  In  general,  equations 
(A4-4)  -  (A4-9)  have  to  be  solved  backwards  in  time.  For  this  purpose  we  make  the 
substitutions: 


Let:  T  =  tf  — t,  -»dT  =  -dt; 


§(t)  =  |(tf-T)=ri(T);  ad  (t)  =  ad  (tf -T)=  ad  (t);  1  =  1,2; 


y  =  x,y,z.  Lienee  the  above  equations  (A4-3)  -  (A4-8)  may  be  written  as: 


Pl4 

rll 

P25 


dru 

dT 

d4  2 
dT 

d4  3 
dT  r33 


9  4  P 14 1 


(ad  -ad  ) 

V  XJ  X2  / 


-»»(«*„ -.‘J 


(A4-9) 

(A4-10) 

(A4-11) 
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dri4 


dT 

dh5 


,  P44  |  d  d 

-m  + - 4  -P44l«xi  -a 

ru  '  xl  x2  - 


- =  _T|  2  +  ^^-T| 5  - p55lad  -  ad 

dT  r22  '  51  y2‘ 


dr) , 


dT 


.  P  66 

_T1  3  + - h  6 

r33 


((ad  -  ad  ) 

;(ad  - ad  ) 
*'  yi  >'2  / 

P«.(az,  -<) 


(A4-12) 

(A4-13) 

(A4-14) 


These  equations  satisfy  the  boundary  condition  that  r\  (d)  =  ^(t ,)  =  0 ,  and  must  be 
solved  backwards  in  time,  that  is,  T  ->•  0 .  We  shall  regard  r|  as  time-to-go  equivalent  of  2, . 


A.4.1.  Analytic  Solution  of  the  VRDE  -  Special  Case 

Analytical  solution  of  the  VRDE  is  possible  for  the  case  when:  s„ 
Si4  =s25  =s36  =s2  =0;s44  =s55  =s66  =2s3  =0  and  rn  =  r22  =  r33  =  2r .  For 
section4.1): 

6rS| 


P11  =  P22  =  P33  - 


Pl4  -P25  -P26 


P44  -  P55  -  P66 


[ar  +  s,  T3] 
6rs,T 


3r +  s 


lT3] 


6rs,Tz 
3r  +  s,  T3 


-  s22  -  s33  -  2sU 

this  case  (see 
(A4-15) 
(A4-16) 

(A4-17) 


Multiplying  both  sides  of  equations  (A4-9)-(A4-ll)  respectively  by 
we  get: 

>44^ldTll  P44 

—  h4  —  P44 


fn  \ 

P44 


VP14  J 


P55 

VP25y 


P  66 
P  36 


P44  ^ 

dhl 

Pl4  , 

dT 

P55  ^ 

dh  2 

P25  j 

dT 

P44 

drl  3 

Pl4  , 

dT 

dtU 

dT 


(ad  -  ad  ) 

=  £fl1l5_P5S(ad  - ad  ) 
Jr  13  33  v  yi  Y2  ' 

(ad  -ad  ) 

V  ZJ  Z2  / 

(ad  -  ad  ) 

V  xi  X2  / 


2r 
P  66 

—  h6-P66v  Z1 


P44 

-hi  +— h4-P44,-xl 


drl3  P55  /  d  d  ^ 

-  =-h2+irh5-P55layi-ay2J 

(ad  -ad  ) 

V  zi  zi  / 


dT 

dh6 


,  P  66 

=  -T1  3  + - h  6  _  P  66 

dT  2r 


(A4-18) 

(A4-19) 

(A4-20) 

(A4-21) 

(A4-22) 

(A4-23) 


Subtracting  equations  (A4-18)  -  (A4-20)  respectively  from  equations  (A4-21)  -  (A4-23) 
and  rearranging  the  terms,  we  get: 

dTl  4 

dT  -1Vl  Ip, 4  I  dT  "  dT  dT 
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dr| , 


dT 
d4  5 

~dT 


-  "H  2  + 

=  43  + 


^P55  1  drl 


VP25y 
r 


ri2  +T^-  =  — (Tti2) 
dT  dT  dTv  ' 


P66^|  dq3 


=  ri  3  +T^-  =  —  (Tti3) 
vp36J  dT  dT  dT'  IJ' 


which  gives  us: 

T1 4  =  Trl  1 ;  ri  5  =  Trl  2 ;  6  =  Trl  3 

Substituting  from  equation  (A4-24)  back  into  equations  (A4-21)  -  (A4-23), 
for  pjj  from  (A4-16)  -  (A4-18)  gives  us: 


dr| 


3s.T 


dT  [sr  +  s,  T3]  [ 3r  +  sx  T3] 

11,12  3S|X!  6rs‘T  (a"  -a") 

V  yi  yi! 


6rsjT 


(ad  -  ad  ) 

\  XI  x2  / 


dT 

d4  3 
dT 


1  2  + 

3r  +  Sj  T3 J  [3r  +  SjT3 


3s3T 


6rsjT 


ll  113 

3r  +  Si  T3J 

3r  +  sj  T3 

(ad  -ad  ) 

\  zi  z2  / 


d  ' 

(3r +  Sj 

(ad 

-ad  )= 

=  (ad 

> 

-ad 

Id 

^rsjT2 

dT 

\  u 

x2  J 

'dT 

A[ 

^3r  +  Sj 

T3)42]=6rSlT( 

ad 

-a"  )= 

=  (ad 

-  ad 

)d 

(3rSlT2 

dT1 

*  yi 

y2 ' 

V  yi 

y2 

'dT 

A[ 

^3r  +  Sj 

T3)r,3]=6rSlT( 

ad 

-  ad  )= 

=  (ad 

d  1 
-  a  j 

|d( 

3rsjT2 

dT1 

k  zi 

z2  > 

\  zi 

z2  > 

'dT' 

Assuming  ax.  ,ay.  ,az.  are  constants;  equations  (A4-28)-(A4-30)  give  us: 
hi  = 

r(2  = 


3rsjT2 

(ad 

V  X! 

-  ad  ) 
x2  > 

3r  +  SjT3 

3rs!T2 

(ad 

V  yi 

-») 

3r  +  SjT3 

3rSiT2 

(ad 
'  zl 

-  ad  ) 
z2  > 

3r +  s,T3 

43  = 

And  it  follows  from  (A4-24)  that: 

44  = 

4s 

46 


3rs  jT3 

(cxd 

V  xj 

1 

P 

x  a 

K) 

3r  +  SjT3 

3rs  jT3 

(ad 

v  yi 

C4 

■a 

s 

i 

3r +  SjT3 

3rsjT3 

(ad 

V  Z! 

-ad  ) 
z2  ' 

3r  +  SjT3 

(A4-24) 

substituting 

(A4-25) 

(A4-26) 

(A4-27) 

(A4-28) 

(A4-29) 

(A4-30) 

(A4-31) 

(A4-32) 

(A4-33) 

(A4-34) 

(A4-35) 

(A4-36) 


UNCLASSIFIED 


37 


DSTO-TR-2838 


UNCLASSIFIED 


Finally,  the  disturbance  term  in  the  feedback  guidance  may  be  written  as: 


Kfn  =  Rr1BTn  =  -1- 

-  2ri 


b  4 
bs 
b  e 


1 

2ri 


arsjT3 
3r  +  2s!T3 


k%  =  R2 


2r-> 


b  4 
b  5 
b  6 


1 

2r2 


SrSjT3 
3r  +  2s!T3 


(ad  -ad  ) 

(ad  -ad  ) 

V  yi  J2> 

(ad  -ad  ) 

(ad  -ad  ) 

(ad  -ad  ) 

(ad  -ad  ) 
\  zi  z2  / 


(A4-37) 


(A4-38) 


A.5  Sight  Line  Rates  for  Small  Angles  and  Rates 


In  order  to  establish  connection  with  between  the  optimal  guidance  and  the  PN  and 
APN  we  shall  assume  that  the  engagement  trajectory  is  such  that  the  azimuth  and 
elevation  sightline  angles  (vp  2 1  ^  ©  2 1 )  (Figure  A-l)  remain  small  during  engagement,  that  is, 
the  trajectory  remains  close  to  collision  close  geometry.  For  this  condition  it  follows  that 
the  interceptor/ target  relative  velocity  is  pointed  approximately  along  the  sight  line  and  is 
approximately  equal  to  the  closing  velocity  Vc .  Under  these  conditions: 

xi2  «r12  »VCT;  ^-x12  =u12  «— Vc  (A5-1) 

at 


Where: 
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ri2  =\xi22  +yi22  +zi22)2  ;is  the  separation  range  between  the  interceptor  and  the  target. 
Interceptor  is  designated  the  subscript  (1)  and  the  target  is  designated  the  subscript  (2). 

Also: 

?i2  »  ri2  *  VCT ;  —  r12  »  -Vc  (A5-2) 

dt 

Where: 

Fj2  =  ( x122  +  y  12 2 ) 2  :  is  the  projection  of  separation  range  on  to  the  x-y  plane. 


We  shall  first  derive  expressions  for  sightline  rates  of  the  interceptor  (1)  w.r.t  the  target 
(2)  and  then  proceed  to  convert  these  to  sightline  rates  of  the  target  w.r.t  the  interceptor; 
where; 

(\|/ 12  ?  0i2 )  are  respectively  the  azimuth  and  elevation  sightline  angles  of  the  interceptor  (1) 
w.r.t  the  target  (2). 

(\|/21  =-y2i,  0i2  =-02i )  are  respectively  the  azimuth  and  elevation  sightline  angles  of 
the  target  (2)  w.r.t  interceptor  (1). 


( . 

A 

( 

A 

.  yi2  ■  2 

tan\|/12  = - ;-*\|/12sec  \|/12  = 

yi2 

ynxn 

;^¥i2  = 

yn  xi2 

yi2  Xi2 

2 

x12 

^12 

x12  J 

l  r12 

r12  J 

¥12 


W2  +.  yi2 


VCT  VcT:  , 


¥21  =-¥12 


x12_  +  _  y12 


VCT  VT2 


(A5-3) 


tan0 


(  . 

j_  A 

(  •  - 

-  A 

12  =4^;  ^0i2  sec2012  = 

z12 

z12  r12 

;  4 

012  - 

z12  r12 

2 

z12  r12 

2 

r12 

lF12 

r12  J 

l  r12 

r12  J 

012  - 


Wl2  +  z12 


VCT  y  x: 

V  c  Tc 1  j 


021  —  —  017 


1 2 


Wl2  ,  z12 


VCT 


yj 


(A5-4) 


It  also  follows  from  equations  (A5-1)  that  the  expression: 


VCT 


1  1 

2  x12  +^YU12  =-■ 


(A5-5) 
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